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Dieter Lasser 
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Abstract: We present a B-spline-Bezier representation of r splines, curvature and 
torsion continuous quintics which have been introduced in CAGD by Hagen in 
1985. Explicit formulas are given for the conversion of the B-spline-Bezier represen- 
tation to the r spline representation and vice versa, and conditions and certain ranges 
of tension values are derived which insure the positivity of the design parameters. 



0. Introduction 



In 1974 Nielson [Nielson 74] gave a piecewise polynomial alternative to splines under tension 
[ Schweikert 66], [ Cline 74], the so-called v spline. The v splines are curvature continuous in- 
terpolating cubics and they are the solution of the minimization of 

t /V 

fVwii 2 ^ + Jhjx'rf , v 7 >o (i) 

/= 0 



over the space 



H 2 = {X: X" e L 2 \_Iq, t N ], X' absolutely continuous on [f 0 , t N ] } 



subject to interpolation conditions and certain end conditions. [ Nielson 74. 5(5] 

In [ Boehm 55] a B-spline-Bezier representation of v splines was given, the y spline, and Boehm 
also pointed out the close relation to Barskys /? splines. [ Barsky 5/] Some details on this relation 
can be found in [ Fritsch 86] especially conversion equations between /? splines and v splines. 
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In 1985 Hagen [ Hagen 55] generalized Nielsons approach by using 

fV'W dt + v / > lII x<L W* . v /x > 0 

\ 7=0 L= 1 

(K ;> 2), rather than (1), minimizing now with respect to 

H K = {X: e L 2 [7 0 , t N ], X ^“ !) absolutely continuous on [7q*0v]} 

and satisfying interpolation and generalized end conditions. Hagens' s concept of 'geometric spline 
curves' includes for K = 2 Nielsons v splines and yields for K = 3 to curvature and torsion con- 
tinuous quintics, the t splines. 

The aim of this paper is to give a B- spline -Bezier representation of r splines. We also derive the 
conversion equations between the Bezier and the originally given Hermite representation in \Hagen 
55]. We discuss positivity conditions on the design parameters of the Bezier representation and 
value ranges for the point weights, the v/s. 

Because we like to find a B-spline-Bezier representation of t splines, we first introduce the Bezier 
representation of segmented curves. In section II we give a short discussion of Nielsons v splines 
w'hile Hagens' s t splines are discussed in section III. 



/. Bezier representation of segmented curves 

Let X(7) be a planar or spatial parametrized curve defined with respect to a partition of the domain 
space by 'knots' 



l o < t\ < - t N . 

The parameter space segmentation induces a curve segmentation in Segments 
X/ : [ tj , 7;^j] — *■ ( d = 2,3). A local parameter u € [0, 1] can be introduced such that 

/ = 0, ... , N — 1 X(0 = X/u) for t e [f /f t M ] 



by the linear interpolation of tj and h+v 

t = (1 -u)t : 4- ut I+x , 

The derivatives have to be calculated now by the chain rule, i.e. 

xW « s - # ir x -«“> - 



where u e [0, 1]. 



1 d r 



A/ du 



7*/(“) 



where A , = t M — t, . 



Now the segments might be given in Bezier representation, that means 

n 

X/(«) = £ W^>) ( 3 ) 

k = 0 

where b n/+A e lR rf ( d= 2,3) , u e [0,1] and 

- fljyo -*)"-* 

are the (ordinary) Bernstein polynomials of degree n in u. The coefficients by are called Bezier 
points. They form in their natural ordering given by their subscripts the vertices of the so called 
Bezier polygon. 
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The derivatives of X,(w) with respect to u are given by 



where 



d r 



du 



T X/(“) = 



nl 



(n-r) 



n—r 

Y ^ b m+k B k~ r ( u ) 



k = 0 



A\ = £T'(b a +\-b a ) 



so that in X(f ; ) = X 7 _j( 1) = X ; (0) , the common boundary point of X /M (u) and of Xj(u) , we have 
for the left sided derivatives of X(r) 



X (r) (f; _ ) = lim X (r) (0 
‘7 

and for the right sided derivatives 

X (r) (f/ + ) = lim X (r \t) 




nl 

( n-r)l 



A r b 



nl-r 




n ! 

(n-r)\ 



A r b 



nl * 



II. Nu-splines 

Nielson s v splines are solutions of the minimization of (1) over the space H 2 subject to the in- 
terpolation conditions 

/ = 0, X(tj) = X 7 
and one of the following end conditions 

i. ) x'(t 0 ) = Xq, x'(ttf) = x;, 

ii. ) X'Vo) = v 0 X'('o) , X"(V) = v^X'(f N ), 

iii. ) X(r 0 ) = , X'(f 0 ) = X'([ N ) , X"(t 0 + ) - X'(f^) = (v 0 + v„) X'(«o) • 



v splines fulfill at any knot (I = 1, ... , N — 1) the continuity conditions 



>< 

11 

j* 


(4) 


X’Ut) = X'(fD 


(5) 


x'(f/ + ) = x*(0 + v ,xvD 


(6) 



what can be written in matrix form as 

x; - a,x; 

The (r+ l) 2 matrix A f is called connection matrix and the 'vector Xf with 1 elements of 
R* (d> r, here is r = 2) is sometimes called the r-jet of X. . 



Because the curvature of a planar resp. spatial curve is given by 

|X',X"| 

k = — resp. k = 

m 3 

we see that a v spline is curvature continuous. 



HX' x XU 
IIX'H 3 



(7) 
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A Bezier representation of v splines can be derived by inserting (3) for n = 3 into (4) to (6). We 
get as continuity conditions of the Bezier representation (Figure 1) 



(1 + q,)b, - <?/b 2i-i + b 3/+1 



where 



(1 + V/*?/) b 3/ _j 
(V/ + <7/)b 3/+l 



V/^/ b 3/-2 + S I 
h s i + yi b 3i+2 



and 



4/ 



V/ 




1 + 



1_ 

_J 

1 + 




( 3 ) 



/ = 1, ... , N — 1 . (8) allows the evaluation of the y 7 's of the Bezier representation of a v spline, 
i.e. the evaluation of the y 7 's for given v ; values. On the other side, the corresponding v spline 
to a given curvature continuous cubic Bezier spline, a so-called y spline, has v 7 values given by 

1=1,..., n- 1 v 7 = + (9) 



(9) was first given by Boehm in [Boehm 55]. He also presented a B-spline representation for cur- 
vature continuous cubics, and pointed out the relation to Barsky's uniformly-shaped /? splines. 
[ Barsky 5/] The connection to general /? splines, sometimes called explicit, discrete or discretely- 
shaped /? splines [Hoellig 55], [ Bartels et ai 57], was also pointed out by Nielson and given by 
Fritsch (see \_Fritsch 56]). 1 In fact, looking at the connection matrices of y splines, v splines and 
(general) /? splines \_Dyn et ai 55] we see that (3 , y and v splines are nothing else than different 
representations of curvature continuous cubics. 



s i 




Dyn and Micchelli \_Dyn et ai 55] have shown that the existence of non-negative local support basis 
functions which sum to one follows - for geometrically continuous spline curves, i.e. spline curves 
with continuous differential geometric invariants like curvature, torsion, etc. - from the total 
positivity of the connection matrix. For v splines the total positivity of A ; has the meaning 
v 7 > 0 and that is exactly the range for the tension parameters v ; covered by the minimum norm 



1 The first local basis for GC 2 splines was developed by Nielson and Lewis in 1975. [Lewis 75] A local 
basis for uniformly-shaped p splines was given in [Barsky 5/] and for discretely-shaped ft splines in 
[Bartels et ai 54], see also [Cohen 57]. 
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characterization of v splines. [ Nielson 74, 86] For v ; >0 we obtain from (9) the y ; range: 
y ; ^ 1 (remark: v ; = 0 y 7 = 1 is the usual C 2 cubic spline). But working with y splines we know 
that y ; > 1 are possible as well. Indeed, if we request positive design parameters, i.e. y ; > 0, so that 
such important properties like the convex hull and the variation diminishing property are given, (8) 
yields to 

/ = i, .... — i v i > + <?/)• (10) 

Hence not only positive tension values but v 7 values in the range given by (10) guarantee positive 
design parameters (remark: v 7 < 0 yields to y ; > 1) and therefore properties like the two mentioned 
above. This result goes conform with work done by Barsky [Bar sky 84] who extended the theory 
of v splines by identifying certain ranges for the v/s that guarantee a unique solution of the in- 
terpolation problem. In the special case of a uniform, an equidistant parametrization, i.e. 1, 
Barsky gives the range v 7 > —4 which is also given by (10). 



III. Tau- splines 



Hagens' s t splines are solutions of the minimization of (2) over the space H K for K = 3 subject 
to the interpolation conditions 

/ = 0, ... , A/ X(r 7 ) = X 7 



and one of the following end conditions (L = K , ... , 2{K — 1) ) 

i.) X aK -'- L \i 0 ) - xj ?*' 1 " 0 . , 

li.) X ( ^(fo") = v 0i2 K_i_ L X< 2 * 1 , X (L> (f/v) = v ,\’,2K-\-L^ 2K 1 L \ 1 n) ’ 

iii.) X(<b) = \(t N ) , X (2 *-‘- L) (r 0 ) = X (2K -'- L) (? tV ) , 

X^ L \iq) — X (L) (^) = ( v o, 2 A'-l-L — v N, 2 K -\-0 X 2K * L \ l o) ■ 



t splines fulfill at any knot (/ = 1, ... , N — 1) the continuity conditions 



X(t/ + ) 


= X(fD 


(ii) 


X'('/ + ) 


= X'(t;-) 


(12) 


X"(t, + ) 


= *VT) 


(13) 


x'V) 


= X"'(t/ ) + V; j X'(t/ ) 


(14) 


x'V) 


fZ, 

>c 

> 

1 

i rz 

>c 

ii 


(15) 



Because the curvature of a spatial curve is given by (7) and the torsion by 



| X', X". X'” 1 
HX'xXf 

we see that a t spline is curvature and torsion continuous. 
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The connection matrix of a general torsion continuous contact of two curve segments of 
U d (d> 3) of degree n (n > 4) - we speak about geometric C 3 continuity , briefly GC 3 continuity 
- is given by [Dyn et al. 55], \_Lasser,Eck 55]: 



A 



GC 3 



10 0 0 

0 a)|| 0^ 0 

0 o 12 cofi 0 3 

0 0)23 w \3 “Ml 



Comparing A GC 3 with (11) to (14), we see that r splines do not take advantage of all shape pa- 
rameters offered by the concept of torsion continuity. They rather form a subset of the set of GO 
continuous curves. 

Furthermore we like'to mention that r splines and visual C 3 continuous curves, briefly VC 3 con- 
tinuous curves, i.e. curves having contact of order r [ Geise 62] with r = 3 , are spanning two ' almost 
totally separated subsets of the set of GO continuous curves. The connection matrix of two curve 
segments having contact of order 3 is given by [Dyn et al. 55], [Lasser 55]: 

1 0 0 0 ~ 

0 0 0 

0 l> 2 °1 Qj 

0 U3 3uj02 °i 




and therefore, comparing A^ C 3 with A GC 3, VO continuous curves form a subset of the set of 
GO continuous curves, and, comparing A^ C 3 with (1 1) to (14), we see that the only curves beeing 
r splines as well as VO continuous curves are the usual O continuous curves. 



To fmd a Bezier representation of r splines we insert (3) for n = 5 into (11) to (15) and get as 
continuity conditions of the Bezier representation (Figure 2) 



(1 + qj) bj 

(1 + yfli) b 5 /— 1 
(y l + <7/) b 5/-H 

(1 -f Stf/) b 5/ _2 
(6/ 4- tfii) s / 
( e / + ?/) b 5/+2 

U + P/<?/) b 5/-3 
(Pi + °l c ll) Q l 
(o, + T t q,) q 

( r I + <?/) b 5/+3 

B-Spline-Bezier Representation of Tau-Splines 



^5/— 1 + b 5/+l 


(16) 


V/^7 b 5/-2 + s 7 


(17.1) 


q,S[ + y, b 5/+2 


(17.2) 


*5/^/ b 5/— 3 + e 7 


(13.1) 


e /<7/ e 7 + 


(18.2) 


<?/ e * + e / b 57+3 


(13.3) 


P/<7/ b 5/-4 + *7 


(19.1) 


*/?/<? + Pi { i 


(19.2) 


T,q,t I 4- ofi 


(19.3) 


Qltf + T / b 5/+4 


(19.4) 
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where 



and 

and 



and 



with 



<7/ = 



A 



/-i 



V/ 



= 1 



<5/ = 



1 + 



1 A l 

r~T~ v /,2 

4 - Hr V ^ 



(It?,) 



£ / = 



1 + 



<// A/ 



( I + <//) 



2 3 



v /,2 



Pi 

°t 



= <7/ 



3tf/ - 7, 

2 1 



T,~ 3 + (7V -5;)-^- 



r / = <7/ 



3 - 7, 



3 <7, — 1 + 



1 



7/ = 



( 1 + (?/) 



3 24 7 -' 



1 + 



( I + <7/) 



A / 

y — v /.2 



s, = 



r , = 



, v . 1 - <7/ A/ 

2(i-<?/)+ yrr'yi 

(i^-<7/r - 4 

. , 2q ‘ A < 

i + y~r v /.2 

( 1 + <7/7 3 

-) <7/ A/ 3 

3 - <7/ + T - TT V /,1 

0+<7/) 3 ~ 4 

. , q ‘ A > 



( 1 + <7/) 



3 3 V /.2 



( 20 ) 



( 21 ) 



(23) 

(24) 



(20, to (24) allows the evaluation of the design parameters of the Bezier representation of a r spline 
and by this the Bezier representation is given. The design parameters are of course not independent 
of each other any longer. This is obvious because two shape parameters are given by (11) to (15) 
but live are given by (16) to (19). The dependences are 



<5/ 



< 77/ 

1 - (1 -q,)z, 



resp. 






'll 

<7/ + (1 -<7/)<5/ 



(25) 
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*1 




and 



Pi 



1 (1 + 4/3/)(£/ + <5/)cr/ 

- — — ; — resp. or 

°l 3(1 + <7/)cr/ — (d/ + q /«/) 



(<3/ -f qjZj)Sjp j 

3(1 + qjjdjpj — (1 + 4/3 /)( £ /+ 3/) 



(26) 



a i 



( 3 / + 4 / £ /) £ / T / 

3(1 + 4/) £ / T / “ ( £ / + 4/)( £ / + 3 7 ) 



resp. 



T / 



1 ( £ / + 4 /)( £ / + 3 /)< 7 / 

z l 3(1 4 - qi)e I - ( 3 / + 4 /£/) 



(27) 



(£/ -f 4/)3/P/ ( 1 + 4/3 )s/T/ 

T ' - ,(1 + .„»,)«, " - w (2S) 

That means if an interpolating r spline is given, i.e. v M and v I2 values, then the design parameters 
of the Bezier representation of the r spline have to be determined in such a way, that they fulfill the 
equations (25) to (28) And the equations (25) to (28) are valid if the design parameter are deter- 
mined by (20) to (24). On the other side, the design parameters of a torsion continuous quintic 
Bezier spline curve have to fulfill the conditions (25) to (28) to make the quintic Bezier spline beeing 
an interpolating t spline in Bezier representation. That means either 3 ; or £ 7 can be chosen as 
independent design parameter and, let's say in case of 3 ; , zj has to be determined by (25) as t 
£ ; = £/(3 ; ) and in addition either p 7 or or t 7 can be chosen as independent design parameter 
and, e.g. in case of p 7 , and r 7 have to be determined by (26) and (28) as <y ; = <7 7 (p 7 ) and 
r i = T ; (p 7 ). If the design parameters of the Bezier representation are determined in this way, then 
the Bezier spline is equivalent to an interpolating r spline having the property of niiiiimizing (2) 
and we can calculate the v 71 and v L2 that means the point weights defining the jumps of the third 
and fourth derivatives in the knots of the r spline by 



V I, 2 = T-(l + ?/Hy-l) 



resp. 



v 1,2 = 



3 (!+<?/} , 1 



4/ 
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and 



resp. 



v i,\ 



( 1 + ?/) 2 [( 1 + <?/)( 1 - 3 q,) - ( jj + <■]{)( -jj - 3t/ / )J 



'/.l 



24 ( 1 + </ /) 
A; 3 * 






resp. 



v I t \ 



24 (1 +‘h) r, , w ,, ,, . <7/ w <7/ , x "l 

[0+<7/)(<7/-3) - (1+ — K — -3) J 



and so on. 



The r splines form a subset of the set of GO continuous quintics. The B-spline-Bezier represen- 
tation of Tau splines is therefore identical with the B-spline-Bezier representation of GO contin- 
uous quintics which was given in [ Eck 57], [. Lasser,Eck 55], with the restriction that for a z spline 
the design parameters y Jt S /y c /t p h a I and r ; can not be chosen independently of each other as in 
case of the GO continuous quintic spline curve. We rather have to set y ; = 1 and have to choose 
S : , tj , pj, o; and r j according to the dependences given by (25) to (28). 



Let's discuss now the positivity of the design parameters for r splines. The dependences (25) to 
(28) imply the following: 

Because of (25), positivity of = <5 ; (e ; ) needs 



0 < e < 



1 

1 ~ c h 



in case of 
in case of 



£ > 0 

On the other side, positivity of e ; = e ; ((5 7 ) needs 

6 > 0 in case of 



0 < 5 < • 



- <7; 

1 ” <h 



in case of 



qi< l <=> Aj < Aj_ j 

qj > 1 A ; > A ; _, 

qi<\ <=> A ; < A ; _! 
q j ^ 1 v — ' A j ^ A 7 _ j 



Because of (27) and (28), positivity of p ; = p ; (r ; ) and of <?, = <7 ; (t ; ) needs 

( £ / + ?/)(£/ + <5;) 

T > 3(1 4- <7/)S/ • 

On the other side, because of (26) and (27), positivity of p, = p,(o,) and of r, = r^a,) needs 



<5/ + t ijcj~ 

3( 1 + t//) 

and because of (26) and (2S), positivity of a, = o,(p,) and of r, = r t (p t ) needs 



(1 + <?/■)/)(£/ 4- <5/) 
3( 1 + q/)6, 



B-Splinc-Bezier Representation of Tau-Splines 



9 



The minimum norm characterization of r splines works for non- negative tension parameter v M 
and v /2 only. As for v splines we can extend the theory of t splines by requesting the positivity 
of the design parameters. 

Because of (20) and (21), for £,>0 and <5 ; > 0 the tension parameter v /2 has to be within the 
range 



that means 






(l + ft)' 3 
V/ - 2 <7/ A, 



v /,2 > - 



12 



r ,2 > -d + ^/) 2 T- 



if q,> I <=> A, > A ; _! 
if < 7; = 1 « A I = A 7 _, 
'/ <7/< 1 <=> A; < A;_, 



Thus not only non-negative but also certain negative v /<2 values are allowed. 



(29) 



Because of (22) to (24), for pj> 0, <r ; > 0 and r, > 0 the tension parameter v 7l has to be within 
the ranges given by (30). 



^ 7,1 < “ £ 3"[( 1 + ^/) 3 + ^ Al v I 7 l ] 

A/ 



v /,i 



> — 



24 (l + < 7 /) 3 + 2(1 — < 7 / -h <7/) v />2 



2 + 



1 + <fo 3 7,2 



(30.1) 



(30.2) 



First v j 2 has to be chosen such that e 7 > 0 and (5 ; > 0 is fulfilled, i.e. v ;<2 has to be chosen within 
the range given by (29), than v M can be chosen within the range given above. 

For q t = 1 (30) yields for example to 

< v /,i < (8 + A 7 v /2 ) 

A / A/ 



j Bar sky [ Barsky 84] extended the theory of v splines by identifying certain ranges for the v/s that 
guarantee a unique solution of the interpolation problem. This idea allows especially the consider- 
ation of different end conditions. The same can be done for r splines, and is indeed the topic of 
actual research. 

Furthermore the idea of [Salkauskas 75] and [Foley 86,87] of introducing interval weights can be 
picked up to create interval weighted geometric spline curves 2 minimizing 

YjPll*' + ££ v />z JX ( %)|| 2 dt. 

7 = U •**/- 1 7=0 L=\ 

For A' = 3 we get interval weighted r splines. Actually we are working on this too. 

2 These interval weighted geometric spline curves are in general not curvature, torsion, etc. continuous. 
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